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ABSTRACT 

For the last several years, Glenn-HT, a three-dimensional (3D) Computational Fluid Dynamics (CFD) 
computer code for the analysis of gas turbine flow and convective heat transfer has been evolving at the 
NASA Glenn Research Center. The code is unique in the ability to give a highly detailed representation 
of the flow field very close to solid surfaces in order to get accurate representation of fluid heat transfer 
and viscous shear stresses. The code has been validated and used extensively for both internal cooling 
passage flow and for hot gas path flows, including detailed film cooling calculations and complex tip 
clearance gap flow and heat transfer. In its current form, this code has a multiblock grid capability and 
has been validated for a number of turbine configurations. The code has been developed and used 
primarily as a research tool, but it can be useful for detailed design analysis. In this paper, the code is 
described and examples of its validation and use for complex flow calculations are presented, 
emphasizing the applicability to turbomachinery for space launch vehicle propulsion systems. 

INTRODUCTION 



The flow environment in the turbine section of a gas 
turbine engine is extremely harsh and complex, Figure 1 . 

The hot gas flow entering the turbine from the combustor 
is highly turbulent, and in most cases, at a temperature 
above the melting point of the metallic materials in the 
turbine. The very survival of the engine requires that the 
turbine be adequately cooled over the entire range of 
system operation. However, system weight and 
efficiency requirements dictate that the cooling 
configuration selected have a minimal impact on those 
system parameters. Therefore, the turbine designer 
must have an accurate knowledge of the detailed 
distribution of critical flow parameters such as pressure, 
temperature, and heat transfer coefficient, both in the hot pjg 1 q-ypical Turbine Flow Problems 
gas stream and in the turbine cooling flow path. A lot of 
this knowledge is derived from the company experience 

base, but more and more, the fine details come from a computational fluid dynamics (CFD) analysis of the 
design. A turbine engine for a space launch vehicle propulsion system poses a particular problem for the 
designer in that the vehicle is continually accelerating, requiring that the turbomachinery operate 
effectively over a broad range of flow conditions. This adds to the attractiveness of accurate CFD tools to 
extrapolate beyond the current experience base. 


As CFD researchers developed the capability to analyze the complex flow fields in gas turbine 
engine rotating components, the emphasis was always on modeling the aerodynamics correctly, 
providing good representation of the pressure distribution and flow patterns in the machine. However, the 
capability to accurately model near-wall viscous phenomena such as convective heat transfer and wall 
shear stress was lacking, primarily because of the fine resolution of computational grids needed to 
accomplish this. Near-wall phenomena require accurate computations of the gradients of physical 
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variables near the wall, and this requires that the grid near the wall be more closely spaced and more 
nearly orthogonal than required for mainstream variables. 

Since the mid 1970’s, the NASA Lewis Research Center (LeRC), now the NASA Glenn Research 
Center (GRC), has had a very strong CFD research and development effort in place, but it was not until 
the late 1980’s that researchers finally addressed the near-wall issues. The seminal publication in this 
area was presented by Robert J. Boyle, of NASA LeRC, at the 35 th ASME International Gas Turbine and 
Aeroengine Congress and Exposition, Brussels, Belgium, June 11-14, 1990 (1) . Boyle used the quasi 3D 
viscous flow code developed by Chima* 2 ’ as the computational vehicle to test his near wall modeling ideas 
for convective heat transfer and wall shear stress analysis. He successfully validated the model against a 
number of turbine stator and rotor experimental databases. 

However, the high flow turning in turbine flows leads to significant problems with grid quality in the 
near wall region. One of the solutions to this problem came from the work of Professor Andrea Arnone, of 
the University of Florence, and colleagues at NASA LeRC* 3 ’. Arnone’s approach was to construct a grid 
that did not require periodicity in the region downstream of the blade trailing edge, thus allowing a “better” 
grid in the blade passage, near the walls. He called the two versions of this code TRAF2D and TRAF3D. 
Building on this, Dr. Ali A. Ameri and Arnone* 4,5 ’ collaborated to adapt Boyle’s heat transfer modeling to 
the TRAF3D code. 

In the meantime, Dr. Vijay K. Gam built on the work of Boyle and adapted the Chima code for the 
analysis of film-cooled turbine blades . This was followed by a switch to the TRAF3D code to take 
advantage of the improved grid near the walls (7) . Garg’s initial approach for film cooling was to 
concentrate surface grid points over the exits of the film cooling holes and specify a distribution of in-flow 
boundary conditions at those points. This worked well if the number of holes was small, but the penalty in 
additional grid points throughout the computational domain was significant. A new approach to the 
gridding of complex problems was needed. 

This issue was addressed by Dr. Erlendur Steinthorsson and colleagues* 8 ’. Steinthorsson revised 
the TRAF3D flow solver, developing a new code, TRAF3D.MB, that was specially designed to use 
multiblock grid systems. Under this approach, the computational domain is broken down into a topology 
of unstructured, connected blocks, and a detailed, structured grid is created within each block. A good 
description of this new code and a demonstration of its capability to accurately compute heat transfer in a 
complex flow can be found in the work of Ameri and Steinthorsson* 9 ’. In this report they compared their 
computed heat transfer on the tip and blade surfaces of a rotating turbine blade with experimental data, 
with very good agreement. Steinthorsson et at* 10 ’ then extended the multiblock approach to the analysis 
of internal cooling passage flow, successfully simulating the air flow and convective heat transfer in a 
complex “branched duct” geometry that included many features typically found in turbine blade cooling 
passages. At the same time, Rigby et al* 11 ’ applied the new code to a turbine blade cooling passage 
model consisting of a 1 80° turn in a rectangular passage, further demonstrating the advantage of the 
multiblock grid approach. Rigby et al* 12 ’ then added the complexity of ribs and bleed holes to the internal 
passage flow, successfully using the multiblock method to include a detailed arid and calculation in the 
small holes bleeding coolant flow from the ribbed channel. Garg and Rigby* 13 ^ further extended the 
multiblock analysis to show the importance of the within-hole and near-hole physics in relation to heat 
transfer on a film cooled blade, including in the analysis the coolant plenum, the flow through the film 
cooling holes, and the external hot gas stream. Heidmann et al (14) went another step and considered a 
blade with two independent coolant plenums and multiple film cooling holes on the leading edge, suction 
side, and pressure side, including shaped holes. Their predictions of surface heat transfer are being used 
to guide an experimental effort currently under way at NASA Glenn. 

In 1 999, the NASA Lewis Research Center name changed to NASA Glenn Research Center, and at 
that time, the TRAF3D.MB code was renamed the Glenn-HT code. Most recently, a Fortran90/95 version 
of the Glenn-HT code has been completed, taking advantage of the object-oriented capabilities of 
Fortran90/95 to produce a code that is more flexible and efficient than the previous version. 
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In the rest of this report, some code specifics will be presented and examples of the application of 
the Glenn-HT code to some specific turbine problems will be described, followed by a description of some 
planned improvements to extend the envelope of applicability for the code. 

RESULTS AND DISCUSSION 


DESCRIPTION OF THE CODE 


A description of the code can be found in most of the references already cited. In brief, the Glenn- 
HT code solves the full compressible Reynolds Averaged Navier-Stokes (RANS) Equations in a rotating 
coordinate system, using a multi-stage Runge-Kutta scheme to march in pseudo-time, with multi-grid and 
implicit residual smoothing to accelerate convergence. To handle the complex geometries typically found 
in turbines, the code uses multi-block grid systems. As currently configured the code uses a k-uj 
turbulence model, and no wall functions. A recent paper by Garg and Ameri (15) successfully 
demonstrated the use of two additional turbulence models in the code. Up to now, the code has been 
used as a steady-state analysis tool, but the new version has the capability for time-accurate analysis of 
transient phenomena. One question that needs to be addressed relative to the application to a Turbine 
Based Combined Cycle (TBCC) Accelerator is whether the analysis of the turbine needs to be done in a 
transient mode, or is the transient slow enough that the analysis can be done as a series of quasi-steady 
calculations. 

TURBINE MAINSTREAM & SECONDARY FLOWS 

An excellent benchmark data case for heat 
transfer in a transonic turbine blade cascade has 
been published by Giel et al (16,17) . Figure 2 shows a 
3D representation of their measured heat transfer 
data on the blade and endwall surfaces. Although 
this was a linear cascade, the combination of thick 
endwall inlet boundary layers and high blade 
turning angle produced a highly 3D flow field. Garg 
and Ameri' 15 ’ used the Glenn-HT code to compute 
the same configuration, evaluating three different 
turbulence models against the data. A sample of 
the calculations is shown in figure 3. Incidentally, 
all the experimental aerodynamic and heat transfer 
data generated for this cascade, and details of the 


The hot gas leakage flow through the gap between 
the rotor blade tip and the stationary shroud has a very 
strong influence on losses and durability. Ameri and co- 
workers have published a number of papers* 9,18 " 23 ' where 
they compared Glenn-HT calculations of the flow and 
heat transfer in this region to experimental 
measurements, with excellent results. Figure 4, from a 
Glenn-HT calculation, illustrates the complexity of the 
flow in the tip gap region. This figure uses numerical 
flow visualization to show the path of the various vortices 
generated over the blade tip. The surface heat transfer 
levels are visualized by color, with red being the highest. 


geometry, are available on request on two CDs 
from the NASA Glenn Research Center. 



Fig. 4. Glenn-HT Computed Flow Traces and Heat 
T ransfer for a Squealer-Tip Rotor Blade 
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Figure 5 shows a comparison of tip heat transfer experimental and computational results. This 
experiment was run at the General Electric Corporate Research Center, and consisted of a single-blade 
cascade with contoured side walls and a stagnation line splitter plate. The computational domain had to 





Fig. 5. Computational Grid, and Glenn-HT Computed Blade Tip Heat Transfer Compared to Experimental 
Data from General Electric Corporate Research Center. 


include the entire test section geometry since there 
were no symmetry planes or periodic boundaries. A 
comparison of the measured heat transfer on the blade 
tip in the tip leakage gap and the Glenn-HT calculation 
showed excellent agreement as seen in Figure 5. 

Another tip-clearance flow validation case was 
conducted by Ameri and Rigby (22) where they modeled 
cooled tip experiment run by Kim and 
. Figure 6 is a schematic of a typical cooled 
blade tip, where film-cooling holes on the tip provide a 
blanket of cooling air to protect the tip region. Kim and 
Metzger modeled this configuration, isolating the region 
of the blade pressure side and the adjacent tip region 
with film cooling holes. The Glenn-HT calculations of 
Ameri and Rigby showed excellent agreement with the 
measured film cooling effectiveness. 

TURBINE FILM COOLING FLOW 


an actively 
Metzger' 2 ^ 



Fig. 6. Typical Rotor Blade with Tip Film-Cooling 


A very effective way to protect turbine blades from the hot gas is to blanket the blade surface with a 
film of cooling air. However, for analysis purposes, this adds a significant complexity to the problem. A 
major effort is continuing to validate the Glenn-HT code for the analysis of film-cooled components. 

Figure 7, from Garg and Rigby <13) , compares the experimental results of Camci and Arts (24) with Glenn-HT 
computational results. In this case, periodicity considerations allowed the analysis to be conducted on a 
representative slice of the tested blade, containing three film cooling holes around the leading edge 
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stagnation region. Figure 7(a) and (b) show details of the inviscid grid. Note that the coolant plenum and 
film cooling hole pipes were gridded in detail and included in the calculation. Figure 7(c) shows contours 
of the computed heat transfer coefficient distribution in the vicinity of the cooling holes, and figure 7(d) 
compares the computed and measured span averaged heat transfer coefficient distribution, showing 
excellent agreement. 



Fig. 7 VKI Blade geometry and heat transfer coefficients 


TURBINE INTERNAL COOLING PASSAGE FLOW 


Perhaps the biggest challenge to the Glenn-HT code is to 
analyze the internal cooling passage flows in a turbine blade or 
vane. A cartoon cut-away view of a generic cooled blade is 
shown in figure 8. In general these passages are non-uniform in 
cross-section, contain turbulators to agitate the flow, contain 
sharp bends, have local flow injection or ejection, and in the 
rotor blade, are subject to the body forces generated by rotation. 
Rigby, Steinthorsson, and co-workers' 11,25 ' 7) have concentrated 
their efforts on these problems, validating the code against a 
wide range of experimental configurations. Figure 9 shows a 
comparison of Sherwood Number data from an experiment, Park 
et al , for a square cross section, ribbed, rotating passage with 
a 180° turn and Glenn-HT calculations by Rigby for the same 
configuration. This test used the naphthalene sublimination 
technique to measure mass transfer rather than heat transfer, 
and the Glenn-HT code was modified to solve the mass transfer 



Figure 8, Cut-away view of 
generic cooled turbine blade. 
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equations. Note that the color scale is different between the two cases. The computational grid for this 
problem contained almost 1 .5 million cells, distributed in 64 blocks. The agreement is excellent, 
particularly in the extremely complex flow region in the vicinity of the bend. 

Additional validation is underway for some complex real engine cooling configurations under the 
NASA Ultra Efficient Engine Technology (UEET) Program, but final results are not yet available. 
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Figure 9, Comparison of Measured and Computed Sherwood Number for a 
Square, Rotating, Model Cooling Passage with Ribs and a 180° Turn. 


SUMMARY AND CONCLUSIONS 

The Glenn-HT code, a three-dimensional (3D), Reynolds Averaged Navier-Stokes (RANS) 
Computational Fluid Dynamics (CFD) computer code for the analysis of gas turbine flow and heat transfer 
has been described and examples of its validation and use for complex flow calculations have been 
presented. The code is unique in the ability to give a highly detailed representation of the flow field very 
close to solid surfaces in order to get accurate representation of convective heat transfer and viscous 
shear stresses. An extensive body of literature documenting the code development and validation has 
been published over the last few years, portions of which were highlighted in this paper. Applications 
include: a) the main gas path through a turbine, including the fine details of flow and heat transfer in the 
clearance gap over the blade tip and the complex flow interactions of film cooling flows with the main 
stream flow; b) the complex flow in internal cooling passages, including turbulators on the walls, 
injection/ejection of coolant, sharp turns, and rotation; and c) the simultaneous combined interaction of 
internal and external flows. 
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In its current form, the Glenn-HT code can be used to analyze, over the broad range of operating 
conditions, the turbines designed for application to Turbine Based Combined Cycle (TBCC) systems for 
space launch vehicle propulsion systems. At the same time, code improvements and extensions are 
being developed to expand the envelope of applicability. The ability to analyze unsteady effects is 
included in the latest version of the code, but validation efforts are just beginning. There are also efforts 
underway to adapt the code for conjugate heat transfer calculations, including heat conduction in the solid 
structures in the analysis. Improvements to the turbulence modeling in the code are being actively 
pursued, with the ultimate goal to be able to include a Large Eddy Simulation (LES) capability. Efforts are 
also underway to streamline the grid generation process through automated topology generation. 

The current version of the code can be made available for domestic use on request. 
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